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Parametric Uncertainty Modeling 

Motivation: 

• Robust Control Theory & Tools 

- Required Uncertainty Model Structure: 

— ► Separated P-A Form: 

- Computational Efficiency Depends on 

Dimension of A Block 

— ► Minimal P-A Model Desired: 

• Practical Robust Control Applications 

- P-A Model Difficult to Form for Real Parameter Variations 

- No General Systematic Approach for Minimal P-A Modeling 
— ► Multidimensional Minimal Realization Problem 

=> Problem to be Addressed In this Paper 



Parametric Uncertainty Modeling (cont) 

General Problem Definition: 

Given State Space Model of Uncertain System: 




t = A(p) x + B(p) u 
y = C(p) x + D(p) u 


Any Element of A(p), B(p), C(p), D(p) — ► Explicit Function* of p: 
Uncertain Parameters: p = [ p,, p 2 , ... , p„ ] 

Pi-i. s Pi s Pi... ► Pi = Pi. + 81 = Pi. + *iS, , 1 8, 1 5 1 


Form a P-A Uncertainty Model: 


Jits 


A(5) 


Pll Pl2 
P21 P22 


Jil 


P - Constant Matrices 
A(5) - Uncertain Parameters 
A(5) = dlag( 8,1,, 5,1* . . . . 5 m I_ ) 
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Parametric Uncertainty Modeling (cont) 

General Problem (cont): 

Any Element of A(p), B(p), C(p), D(p) -► Explicit Function* of p 


"Explicit Functional Forms: 
Linear Function" 
Multilinear Function 

Rational Function 


Example : 

»ij(p) = Pi + Pi a. 

**j(P) = Pi + P 1 P 2 «. 

a,.(p) = Pi ± P 2 «. + P^Pi 
IJVF ' Pi Pi + «i p* 


** Formal Solution by Morton & McAfoos (1985 ACC & CDC) 

=> Many Practical Problems: 

Multilinear (Rational, ... ) 


Objective 

Develop : Systematic Method for Obtaining a P-A Model 
Given : State-Space Model of a MIMO Uncertain System 
such, that : 

• Any Element of A(p), B(p), C(p), D(p) Is a Multilinear Function of p: 

aij(p) = pi + pip 2 ao 


• The Resulting P-A Model Is Minimal (or Near Minimal) , l.e.: 

=> A(6) = dlag( 5,1,, Mi, • • • , Kfm ) 

has Minimal Dimension for the Given State-Space Model 

Extend : Multilinear Results to Rational Case 
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General Solution Framework 


Block Diagram Perspective: 


|~u] | A(p)B(p)| fy 1 

'mm 


[J] = S(p)[;] = (S(p„) + S a (6))[*] 


S 4 ( 5 ) 




A 4(8) B a (8) 
C a ( 5 ) D a (8) 

h 

=> 




[;] J 


A(p.) B(p.) 
C(p.)D(p.) 

w 

♦6 


S(P.) 




[2 



S(p») 



= 48) y A 


A(5) = diag( 8 ,I„ . . . , 5„I„ ) 


General Solution Framework (cont) 

Equating Given & Desired Models: 

ills = s (p.) 

A( 8 ) = diag( 8 , 1 ,. 8 2 l 2 , .... S„I m ) 
l|l-A( 8 )f§-«A( 8 )f|f = S a ( 8 ) 


Solution of P21, P12, & Pn Matrices: 


Hf( I - A( 8 ) Hi- 1 A( 8 )i| = S 4 ( 8 ) 

\ f / 

Unknown Matrix Elements 


A 4 (8)B a (8) 
C 4 (8)D a (S) 

t 

Known Matrix Elements 
(Function of 5's) 


General Solution Requires: Direct Matrix Inversion 

( i - a( 6) p„ r l 

=* Symbolic Matrix Inversion & Subsequent Solution 
Difficult for Many Practical Problems 
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Multilinear Solution Framework 


PK I -4(8) H)-> 4(8)11 = 

a f r 

Unknown Matrix Elements 


m = s 4 (S) = 


A 4 (5)B 4 (5)1 


lC 4 (8)D 4 (8)J 

Known Matrix Elamants 
(Multlllnaar Function of 5'a) 


Flnlta Powar Sarias (Exact Solution): 

(I-4(5)P„)-‘ = I + (4(6) P„) + (4(8) P,,) 2 + ... +(4(8)P„) r 

such that: (4(6) P U ) MI = 0 =» Regulres Special Structure for Pn 

where: r - Determined by Maximum Crossterm Order In A, B, C, D 


S 4 (6) = 



= P« [ I + 4(6) P„ + (4(5) P,,) 2 + . . .+ (4(8) P„) r ] 4(8)P l: 


Uncertain Parameter Linear Terms Uncertain Parameter Crosst 

Note : 1.) nth Order Terms 2.) Inverse Terms 

-* Repeated Parameters -> Redefine Parameters 


Ex.: p,* = p,p w 


Ex.: 4- = Pi 

Pi n 


Uncertainty Modeling Procedure 

To Obtain a Minimal (or Near Minimal) 

P*A Uncertainty Model: 

0. Determine Pa and Extract S 4 (5) : 

m - I . - f 4 -® ,, -“ > l 

C( Po ) D( Po ) J [c a (5)D a (6> 

1. (Wins A Matrix: a<5) = diag( SJL l9 ZJL 2 SJ. ) 

Repeated Parameters Ontylot 
nth Order Uncertain Parameters 

2. Determine P 2 t and P 12 Using Linear Terms (Morton & McAfoos): 

mmm= [s 4 (8)) 0 = [ Aa(5)Ba(5) 1 

k C A (8)D A (6)L 


Known Linear Uncertain Parameter Terms Only 
(No Uncertain Parameter Crossterms) 
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A(p) = 


=» A(p) 


Modeling Procedure (cont) 


l Determine Pn Using Uncertain Parameter Crossterms: 


P2 1 (A(5)^ 1 A(5)P 1J = fS 4 (8)J, = 
= [S 4 (8 )]j = 

P al (A(6)||) f A(8)P IJ = [S 4 (8)] p = 


A 4 (8)B 4 (8)1 

C 4 (8)D 4 (8)J, 

A 4 (8)B 4 (8)1 

A 4 (8)B 4 (8)| 

C 4 (8)D 4 (8)L 


Known 
First -Ordsr 
Crosstsrms 


Known 

Sscond-Ordsr 

Crosstsrms 


Known 

rth-Ordsr 

Crosstsrms 


with Nllpotency Condition Satisfied. 

H Pn Cannot be Found such that ALL of the above 
Equations and Condition are Satisfied: 

a. ) Determine which Parameters Need to be Repeated 

b. ) Repeat Procedure from Step 1 Augmenting A Matrix 


Once Pn has been Determined, 

Minimal (or Near Minimal) P-a Model Has Been Found 


Example 


Given Uncertain System Model : 


V 

° ° 


0 

— 0 


0 

WPS it 


C(p) 

•[;::] 


-V, 

m o 

0 

0 

-Vj i 

0 

0 


•vll 


where: Lu 

HJ 

II 
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Example (cont) 

P-A Model Solution: 



Example (cont) 

P-A Model Solution (cont): 

-1 0 0 0 L.. 0 

0 0 -V, V. 0 0 0 

p 2i = 0 0 0 0 1 0 E.. 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 
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Example (cont) 

P-A Model Solution (cont): 



ooo 



0 0 25 ^, 


ooo 
0 0 0 


A(5) = diag [ 5j^ I 2 5^ I 3 8 0 . 8^] 


Extension to Rational Case 
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Extension to Rational Case (cont.) 

System Equations: 



+ 



- Sn.Sp .- 1 



«A 


y A = [Pi2,p.2.]Sd/'[5] + 


Pll “[P12, ^ 12 .] ^ D ,* 1 



Ua 


«A = Ay 4 



Concluding Remarks 

Multilinear Solution Framework 

- Solves Multilinear Parameter Cast 

=> Accomodates nth Order and Inverse Terms 

- Eliminates Symbolic Matrix Inversion In Computation of Pi i 

=> Computationally Tractable for Symbolic Solution 
(Symbolic Algebra Tool Required) 

- Can be Extended to Rational Parameter Case 

=> Preliminary Results 

Systsmstis. Procedure for (Near) Minimal P-A Modeling 

- Minimality Is Relative to Given State Space Realization 

=> A Lower Dimension P-A Model May Exist for Different Realization 

- (Near) Minimality by Construction 

=> Minimality may not Always be Assured 
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Further Work 


• E valuate/Reflne/Generalize Procedure 

- Wider Class of Problems 

- Multidimensional System Theory 

• Automate Modeling Procedure 

- Mathematlca/Maple 

- Output Files to Matlab 

• AddIv to HSCT Problems 

- Configuration Evaluation 

- Control System Analysis & Design 
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Abstract 

Advanced robust control system analysis and design is based 
on the availability of an uncertainty description which separates the 
uncertain system elements from the nominal system. Although this 
modeling structure is relatively straightforward to obtain for multiple 
unstructured uncertainties modeled throughout the system, it is difficult 
to formulate for many problems involving real parameter variations. 
Furthermore, it is difficult to ensure that the uncertainty model is 
formulated such that the dimension of the resulting model is minimal. 
This paper presents a procedure for obtaining an uncertainty model for 
real uncertain parameter problems in which the uncertain parameters 
can be represented in a multilinear form. Furthermore, the procedure is 
formulated such that the resulting uncertainty model is minimal (or near 
minimal) relative to a given state space realization of the system. The 
approach is demonstrated for a multivariable third-order example 
problem having four uncertain parameters. 

1. Introduction 

Advanced robust control system analysis and design is based 
on the availability of an uncertainty description which separates the 
uncertain system elements from the nominal system. More 
specifically, the uncertain system components are contained in a block- 
diagonal A matrix, which is connected to the nominal system, P(s), 
such that the closed-loop uncertain system is described by a linear 
fractional transformation (LFT). The idea of separating the uncertain 
part of a system from its nominal part in this manner, for use in robust 
control system analysis and design, was first posed by John Doyle (see 
[3] and [4]), and the robust control theory associated with this 
structured description of uncertainty continues to be an important area 
of research. A block diagram of this modeling structure can be 
depicted as follows in Figure 1: 



Figure 1. Block Diagram of General Uncertain System 

where u contains all external inputs to the system (e.g., disturbances, 
control inputs, etc.), y contains all outputs from the system (e.g., 
controlled outputs, measured outputs, etc.) and u^ and connect the 
uncertainties represented by A to the nominal system, P(s). Although 
this modeling structure is relatively straightforward to obtain for 
multiple unstructured uncertainties which occur throughout the system, 
it is difficult to formulate for many problems involving real parameter 
variations. Furthermore, it is difficult to ensure that the uncertainty 
model is formulated such that the dimension of the resulting model is 
minimal (i.e., the number of repeated parameters in A is minimized). 
Although formulating an uncertainty model is a requirement for 
utilizing the recently developed robust control analysis and design 
techniques mentioned above, very little research has been reported in 
the literature which addresses this problem, particularly for the real 
parameter uncertainty case. Results to date primarily apply to multiple 
uncertain parameters which enter the system model in a linear 
functional to rm, although some work involving nonlinear special cases 
have been worked [10]. The results for linear uncertain parameters 
were first presented in [8] (Morton & McAfoos, 198:5) and [9] 


(Morton, 1985). A later paper [10] (Steinbuch, et. al. t 1991) 
summarizes the general uncertainty modeling problem and the results to 
date, and presents two simple scalar nonlinear uncertain parameter 
examples. However, no solution to the general minimal uncertainly 
modeling problem has been found. The objective of this paper is lo 
present an important extension to these uncertainty modeling results. 
Specifically, a procedure is presented for obtaining a minimal (or near 
minimal) uncertainty model (having the form of Figure 1) given the 
state space realization of an uncertain system with multiple parametric 
uncertainties entering the model in a multilinear functional form. It 
should be noted that minimality here is relative to the given state space 
realization. As discussed in [1] and [2] (Belcastro, et: al., 1989 and 
1991), the dimension of the uncertainty model (i.e., the dimension of 
the A matrix) is dependent on the state space realization of the system. 
Thus, one can consider the minimality of an uncertainty model for a 
particular state space realization, or one can consider the achievable 
minimality of the uncertainty model irrespective of the system 
realization. In this paper, we present a method of obtaining a minimal 
(or near minimal) uncertainty model relative to the given state space 
model of the uncertain system for multiple uncertainties entering the 
model in a multilinear functional form. The multilinear frameworit 
significantly reduces the computational complexity involved in 
obtaining a solution, as compared to solving the problem directly for 
the rational parameter case. Moreover, it can be shown that the 
multilinear solution framework can actually be used to solve the 
rational parameter case, as well. Thus, it provides a means of 
determining an uncertainty model for many difficult problems of 
practical interest. 

The paper is organized in the following manner. Section 2 
presents a formal problem definition for the general uncertain parameter 
case, briefly summarizes results for the special case of linear parametric 
uncertainty, and defines the problem to be addressed in this paper. 
Section 3 summarizes our results for this defined problem, and Section 
4 presents an example problem which demonstrates these results. 
Section 5 briefly discusses the application of the multilinear solution 
framework to solve the rational uncertain parameter problem, and 
concluding remarks are given in Section 6. 

2. Parametric Uncertainty Modeling: 
Problem Definition 

2.1 General Problem Definition 

Consider the state space model of an uncertain system: 

x = A(p) x + B(p) u , x e R n \ u e R"* (la) 

y = C(p) x + D(p) u , y € R"* (lb) 

where p represents a vector of real uncertain parameters: 

p = (Pi.P2. --.Pml e R m (2) 

It is assumed that each entry of the model presented in equation (1) is a 
function of the parameters p. For the general rational case considered 
in this paper, the uncertain parameters can appear in a rational 
multivariate functional form within each element of the system model. 
For example, as given in [10] (Steinbuch et. al., 1991), the (i j)^ entry 
of the A matrix could have the form: 
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Aij(p) = 


Pi + P2 a<> + P2 2 P3 
Pi P3 + »I P4 


(3) 


where ao and are constants. It should be noted that n^-order terms 
are included here because they can be handled within a multilinear 
framework by defining n-1 additional uncertain parameters which are 
equal to the parameter being raised to the n* power. For this example, 
a new uncertain parameter, p2* ~ pj, could be defined and p2 2 would 
then be replaced by p2p2*- 

f The uncertainty modeling problem consists of three 
components: scaling of the uncertain parameters, extraction of the 
uncertainties from the nominal system, and formulation of a linear 
fractional transformation (LFT) (see [5], Doyle, et. al., 1991 for a 
review of LFTs). These components are reviewed below. 

Uncertainty Scaling: 


y A = p n u A + p n[ *] 

(13) 


(14) 

U A = A(8) y A 

(15) 

A(8) = diag(8il,. 82I2 6ml m ) 

(16a) 

A(S)€ R T, a* n 4 

(16b) 

TT» 

n A = Z r i . *1 = dim© 

(17) 


i-1 


Each uncertain parameter pi in p can be bounded by an upper 
bound, Pmaxi« and a lower bound, pminji as follows: 


Pminj ^ Pi ^ Pmaxj W 

Then the parameter can be written in terms of some nominal value 
within this range of uncertainty. One way to do this is shown below: 


Pnomj 


+ Pnorrij + Sj 8j 

(5) 

_ Pminj ^ Pm«j 
2 

(6) 

Pmax- " Pminj 

2 

(7) 

I8jl £ 1 

(8) 


Equations (4) - (7) can also be written in vector form by stacking each 
associated parameter quantity into vectors. The 8* terms as defined in 
equations (5) and (8) are the uncertain terms that will be separated into 
the A matrix of Figure 1 . 

Uncertainty Extraction: 


where Pn, P| 2 . P 21 . and P22 are constant matrices with P22 = 
S(pnom)' and the matrices Pn. P12. and P21 are related to Sa(S). The 
Ii terms in equation (16a) represent the identity matrix with dimension 
equal to the repeatedness of parameter 5j. For example, the squared 
uncertain parameter of equation (3), i.e. P2 2 , results (after scaling) in 
the term bfi. Thus, this example would require that both $2 a™* $2’ = 
62 (associated with the uncertain parameter pj' discussed above) 

appear in A, which means that I2 in equation (16a) would be a 2- 
dimensional identity matrix. 

The objective of the uncertainty modeling problem is to find the 
matrices Pn, Pi 2. and P21 such that the system of equations 
represented by (13) - (16) is equivalent to the system represented by 
equation (9). To do this, equations (13) - (15) are combined such that 
ua and y^ are eliminated, as follows: 



IP22 + Pzid - A(5)P„) 1 A(5)P 12 ( * ] 


(18) 


Thus, the uncertainty modeling problem can be thought of as a multi- 
dimensional (minimal) realization problem defined by the following 
equation: 

S A (S) = P 2 i(I - A(5)P„X 1 A(6)P,2 (19) 


Using equation (5), the state space model of the uncertain 
system given in (1) can be rewritten in compact form as follows: 


where 8 represents the uncertain parameter vector defined in equation 

( 10 ). 


y ] - S<p) [ * ] = Sftwof * ] + S 4 (8)[ * ] 

where: 

5 = [5|.fe.-.<ta] € R m 


“ C(p) D(p) “ + s a(5) 


S(Pnom) — 


A(p„om) B(p nom) 

C(Pb om ) D(pn om) 


Sa(5) = 


A a (5) 

. Q(8) 


B A (8 ) 
D a ( 8) . 


2.2 Summary of Results for Linear Parametric 
(9) Uncertainties 

As indicated previously in this paper, uncertainty modeling 
results have primarily focused on the special uncertainty case involving 
(]0) multiple uncertain parameters that enter the system model linearly. 
Results for this case were first presented by [8] (Morton & McAfoos, 
1985), and involve solving equation (19) with Pn = 0. For this case, 

P21 and P 12 can easily be found by expanding Sa( 8) as a linear 
combination of the 5j terms, and decomposing the resulting coefficient 
(12a) matrices. If any of the coefficient matrices has rank greater than one, 
then the associated 5j term must be repeated in A a corresponding 
number of times in order to perform the decomposition. For example, 
(12b) if the coefficient matrix for 5j is rank 2, then 5j must appear twice in 
the A matrix. This is also discussed in [9] (Morton. 1985). 


Separation of S(p) into nominal and uncertain parts, S(pn 0m ) and 
Sa(5), respectively, results in the extraction of the uncertainties from 
the nominal system. 

Formulation of a Linear Fractional Transformation fLFTt: 

Equation (9) can be rewritten in the form of an upper (time 
domain) LFT by defining an input vector, ua, and an output vector, 
yA, associated with the uncertain part of the system as follows: 


2.3 Specific Problem Def inition for this Paper: 

Multilinear Par ametric Uncertainties 

In this paper, we consider the case of multiple uncertain 
parameters which enter any element of the system described in equation 
(1) in a multilinear manner. It should be noted that rational multivariate 
elements involving only one denominator term can be represented in a 
multilinear form directly. For example, 
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■ Pl 1 B „V flPi 

P1 ** ( 21 ») 

= P3 + Pl P2P3 *o + Pl P2 2 

*-*•*■* 

redefined. For example, in uncertain model element represented by 
equation (3) could be approximated in a multilinear form as follows. 

a j n \ - Pi + Pi«°- t -P^Pi 

A ^>" PiP3 + ».P4 (20l) 

= ?4 (Pl + P2 *o + P2 2 P3 ) 

WherK ?4 * 1 (20b) 

P4 Pl P3 + ®1 P4 

Thus in this formulation the fourth uncertain parametcrJP* , is 
dependent on the uncertain parameters pi , p3* and p 4 . ms approach 
therefore poses a slight ruction to the general case. Hwmm 
brief discussion of a technique for formulating the rational problem m 
such a way that the multilinear solution framework can be used is 
presented in Section 5. 

2.4 Formal Froblem Si a tefflfiBi 

A formal problem statement based on the above discussion can 
be summarized as follows: 

Given : An uncertain system in state space form as in equation (1), 
i.e.: 

* * A(p) x + B(p)u, x€ R* u€ R"* 
y = C(p) x + D(P) » . y € R % 

which can be rewritten as in equation (9), i.e.: 

[Jl-wIil-sMil * s * , 6 ) li] 

Find : The matrices P2i» Pl2» ^ Pit 5uc h * at * c above sy$ icm can 
be expressed as in equations (13-16), i.e.: 

y A = p n u A + p > 2 [ul 


y ] = P2 ' U A + P “[ u] 


u A = A(8) y A 

A(6) = diag( Sjli, 52b. * • • » Um ) 

A detailed discussion of a solution to this problem for uncertainties 
which are represented within a multilinear framework, as discussed 
above, will be presented in the next section. 


Aa(5) Ba(5) 

Sa( ) = [ C a (5) D A (8) 

= P 21 (I - A(5) p n)-' A(8)P, 2 (22) 

= [ p 2W ] a ' A * 5 > P '»>'' A(S) 1 PlI « P,2 « 1 

where the A a (5), B a (6). C a (8), and D A (5) terms in equation (22) are 
formed by scaling the uncertain parameters p and extracting the 
uncertain 5 terms from the nominal system, as discussed in Section 2, 
and P21 and P12 are partitioned appropriately. Thus, the matrices 
Aa(8), B a (5), C a (6), and Da(6) are known matrix functions of the 5 
parameters, and the matrices P21, P12. p H unknown 
variables for which equation (22) is solved. This section presents the 
main result of the paper - namely a solution to the above FoWanfor 
uncertainties that are represented within the multilinear framework 

described in Section 2.3. . . 

As stated above, the solution to this problem involves solving 
equation (22) for P 2 i. P12. «nd Pji. However, the inversion of the 
quantity ( I - A(8) Pn ) in equation (22) for multiple pwameter 
problems can become very cumbersome because Pn is of the same 
dimension as A(8). and the inversion has to be performed 
symbolically. Moreover, each dement of Pi i must be detomined such 
that equation (22) is satisfied. Within the multilinear framework, 
however, this quantity can be replaced by a finite senes. To see this, 
consider the matrix equation: 

X . A n+ i = (I - A) (I + A + A 2 + A 3 + ... + A n ) 

which can be written for any matrix A. Assuming that the matrix 
(I A ) is invertible, this equation can be rewritten as: 

(I - A)" 1 (I - A"* 1 ) = I + A + A 2 + A 3 + ... + A" 

If matrix A is structured such that A n+1 = 0 (i.e., A is nilpotent) , 
then: 

(I _ A)- 1 = I + A + A 2 + A 3 + . . . + A“ 

This development is similar to the Neuman aeries 

in [6] (Halmos, 1974) for a matrix A such that II A II < 1. For our 

problem, however. A = A(8) Pn. where A(8)isa diagonal matrix and 
Pi i is unknown. Although A(8) is norm -bounded by unity. Pn is not 
norm-bounded. However, since Pn is to be determined, requtrmg p n 
to be structured such that 

(A(8)P U / +l = 0 (23) 

yields; 

a - d(8) Pn) 4 = I + (A(5) Pn) + (A(8) PnX 

+ ... + (A(8)Pn) r (W) 

Substituting this into equation (22) results in: 


Sa(5) = 


A a (8) B a (8) 

C A (S) D a (8) 


= P 2 i[I+A(8)Pii+(A(8)Pn) I +...-KA(8)Pii) r lA(8) p n 


which can be rewritten as: 


As indicated in Section 2, the solution to the uncertainty 
modeling problem posed above involves finding the matrices P2I1 “12* 
and Pn such that the Sa(6) matrices given by (12) and (19) are equal, 


S a ( 5) = P2lA(5)Pi2 + 

P 2 ,[A(8)P,i+(A(8) Pi ,)*+... 

...(A(8)P u ) T ]A(8)Pn 
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The first term on the right side of equation (26) represents the linear 
uncertain components of S A (8), and the second term adds in the 
nonlinear terms. Furthermore, since the nonlinear terms of S^(8) 
consist of cross terms and n*-order terms (which can be represented as 
cross terms), die order, r, of the highest term in the series of equation 
(26) is defined by the highest cross term order required to realize 

S a ( 5). Thus, r is defined by the order of the highest cross-term 
occurring in A a ( 8 ), B a ( 8 ), C a ( 8), and D A (8), i.e.: 

r = max ( O^, Og, Oc, Oj) ) (27a) 

uid Q Af Og, Oc. and Of) represent the order of the highest-order 

cross-product term in A a ( 8 ), B a ( 8 ), C a ( 8). and D A (8), respectively. 
That is, fora general uncertain mxn matrix M: 

Om = max [ order ( ) ; for all i = 

J and j = 1,2,...*] (27b) 

where the order of each my is the order of its highest-order cross- 
product term, and cross-product term order is defined as: 

order ( 8j 82 S3 = i-1 

for i = 1 , 2, . . . , n A . (27c) 

Thus, the maximum value of r is r mix = n A -l, where n A is the 
dimension of die A matrix and is given by equation (17). The nilpotenl 

requirement of equation (23) for (A(8) Pi 1) can be satisfied if the 
elements of Pn . Pij. satisfy the following structure: 

1. ) Pu = 0; i=l,2 n A 

2. ) If pH * 0 , then for 

I = 1, 2 iia and j = 1, 2, . . . , n A : 

a. ) pji = 0 ; 

b. ) Piei j©i = 0 or pi©2j®2 = 0 OT 

• ■ • Pi®(n A -l)4®(n A -i) = 0 

(28) 

where the symbol represents "modulo n A "addition [7] (Horowitz 
and Sahni, 1978) over the set {1, 2, .... n^}, i.e.: 

• a + b if a + bSn A 

aWD ‘| a + b - n A if a + b>n A 

1 , 1 £b <n A 

and nA is the dimension of A (and, hence, Pjj) as defined in equation 
(17). It should be noted that requiring Pji to satisfy the conditions of 
(28) does not impose a restriction in solving the uncertainty modeling 
problem, but rather it is a means of removing unnecessary freedom in 
determining Pn based on the uncertain system being modeled. Thus, 
(28) assists in the process of solving for Pn. 

Using this multilinear framework, P21 and P12 can be found 

using the linear uncertain terms of S A (8), and Pji can be found using 
the nonlinear terms of S A (8) such that the conditions of (28) are 
satisfied. Thus, the procedure presented in [8] (Morton & McAfoos, 
1985) (and briefly described in Section 2.2) for obtaining an 
uncertainty model for multiple linear uncertain parameters can be used 
to obtain P21 and P12, and these matrices can be used in the second 
right-hand term of equation (26) so that Pn can be determined directly 
using equations (26) and (28). Details of the procedure for doing this 
are presented in [1] and [2] (Belcastro, et. al., 1989 and 1991), and an 
example problem is presented in Section 4 which demonstrates these 
results. 

3 J Uncertainty Modeling Procedure 

Obviously, in order to reduce computational complexity in 
robust control system analysis and design, it is desired to obtain an 
uncertainty model of minimal dimension. As discussed in [1] and [2] 


(Belcastro et. al. t 1989 and 1991), the dimension of the uncertainty 
model is dependent on the system state space realization. These papers 
address the problem of obtaining a state space realization of an 
uncertain single-input single-output (SISO) system (given its transfer 
function) such that an uncertainty model of minimal dimension can be 
determined. For practical multivariable applications, however, it is 
usually desired to retain physical relevance to the problem being 
considered in assigning the states of the system, so that a particular 
state space realization may be preferred. Therefore, given a desired 
state space model of an uncertain system, one would like to be able to 
determine a minimal uncertainty model for this particular realization - 
which may or may not be an overall minimal uncertainty model for the 
system. A procedure to obtain a minimal (or near minimal) uncertainty 
model relative to a particular state space realization (based on the 
multilinear framework presented in Section 3.2) is therefore given in 
this section. 

Given a state space realization of an uncertain system whose 
matrix elements are multilinear functions of die uncertain parameters of 
the system, it is desired to obtain an uncertainty model of the form of 
Figure 1, which has a minimal (or near minimal) number of repeated 
parameters in A. This can be done using the following approach: 

1. Define a A matrix of the form of equation (16) which has only 
those repeated uncertain parameters necessary to realize the nth- 
order uncertain terms in die model, as discussed in Section 2.1 . 

2. Follow the procedure given in [8] (Morton & McAfoos, 1985) 
and [9] (Morton, 1985) for the linear uncertain parameter case 
to obtain P21 and P12 using equations (22) and (26). If 
problems with rank occur in defining P21 and P12* go back to 
step 1 and add a repeated parameter to A, as described in 
Section 2.2. 

3. Once P21 and Pi 2 have been obtained, use the nonlinear 
uncertain terms in equations (22) and (26) to obtain Pi 1 such 
that the conditions of (28) and, hence, equation (23) are 
satisfied. If Pn cannot be determined such that all of these 
equations and conditions are satisfied, the dimension of A is not 
large enough. If this occurs, it must be determined which 
parameter must be repeated (based on the specific problem 
encountered in trying to satisfy the above equations), and the 
process begins again at step 1 with the repeated parameter being 
added to the A matrix. Once Pi 1 has been successfully 
determined such that all equations and conditions are satisfied, 
the minimal (or near minimal) uncertainty model for the given 
state space realization of the system has been determined, and 
equations (13) - (16) can be used to model the uncertain system 
as depicted in Figure 1 . 

It should be noted that the above procedure yields a minimal (or near 
minimal) uncertainty model by construction , since the initial A matrix 
defined in step 1 is of the smallest possible dimension required to 
model the given system, and additional parameters are added to this A 
matrix in steps 2 and 3 only if required. An example problem 
illustrating the above procedure is presented in Section 4. 

4. Example 

Consider the third-order multivariable system described in state 
space form as in equation (1) by the following realization: 


(29a) 
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B(p) = 




0 

0 




(29b) 


A(po) = 


V.U. 

0 

0 


-VJ^ 


-» 1 
CTwoU. 


0 

0 

-v£* 


(35a) 
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where the uncertain parameters L«, Ur, 0u. and O w vary over the 
following ranges: 


B(po) = 




0 

0 


105.7 $ Lo 

£ 841.1 

(30.) 

C(p„) = 

[ 1 0 0 

10.4 S U 

5 795.5 

(30b) 

L 0 0 1 

5.74 £ Ou 

£ 9.69 

(30c) 


r 

3.95 £ a w 

5 13.4 

(30d) 


I -v # Su 


Hie element, of equation (29) can be exprescd as multilinear functions 
of the uncertain parameter, a, follow,: 


A(p) = 


-V.U 

0 

0 


-V.U 
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0 

0 

V.L. 


OwUr 

<^yy 

B(p)=| 0 
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0 0 
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B a (5) = 




0 

0 


-V.(l -^)5u 


where: 


where: _ 

^ = ^ = u 

.001189 S U £ *009461, .0013 £ U S 00962 

The first step is to extract the uncertain 5 terms from the nominal 
system by scaling the uncertain parameters as in equation (5), as 
follows: _ 

U = U* + = U* + 

Ur = Ur # *u«U s 

<J» = <fc. + *o,V. = On. + 8s. 

Ow = Ow. + *oA. = Ow. + 6o. 
so that, as in equation (12): 


(34) 


where: 


2n 

c .<«> ■ i ns ] • d * < »-[;si 

a Al = 2o 2L. t ^u^* 

- 2 T T_ 2 

+ o w .8l. +So.6l. 

b Al — 20 b.Lu.5l. ^ Lb.6q. + 2t^.6L,5o. 

+ o».8l. + $o.5l. 


(36b) 


(36c) 


(37s) 


(37b) 


— 2Cw.Lw.5l. L..&5. + 2Lw.6l.5b. 

+ Ow.8u 2 + $b.£. 2 ( 37c ) 

As can be seen by the last term in equation (37) (for either a^,, b A ,. or 
b^), r = 2 for this example problem (as defined by equation (27)). 

Since S A (S) contains 2 nd -order terms associated with U, and Ur. Ow 
terms associated with these variables will have to appear twice in A 
Thus, the dimension of A going into Step 1 of Section 3.2 is six. For a 
six-dimensional A, the matrices P21 and P12 can be determined, as 
described in Step 2 of Section 3.2. However, it is impossible to obtain 
a Pi 1 matrix which satisfies all of the equations discussed in Step 3 of 
Section 3.2. Moreover, it is determined in that step that the 6 term 
associated with L w must be repeated a third time. Therefore, when 
steps 1 - 3 of Section 3.2 are repeated, the resulting uncertainty model 
can be expressed as in equations (13) - (16) and (22), where: 


A = diag(8u 5, 8u 5. K. 

= diag( Sl.U Su*3 So* * 0 * 1 


P2U = 


-1 aJU. 0 0 0 0 

0 0 -V, V, 0 0 0 

0 0 0 0 1 0 Lw. 


(38) 
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Pl2* = 


Pl2» = 


0 

0 
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0 

0 

0 
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0 

0 


o -fu v * 

0 0 ' o 

0 o 
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S st 

2S Ur® w <,f w <V^ 




Pll = 


0 

0 

L 0 


0 p„ 0 0 0 0 0 

0 0 0 0 0 P M 0 

0 0 0 0 P 37 

0 0 0 0 P 47 

p 53 0 0 0 p J7 

0 0 0 0 0 

0 0 0 0 


0 0 
0 0 


0 

0 

0 
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Pl2= -SEulh 


P26 


.Si 


Pj7= 

P4?= 

P5J = 

P 57 = 2S u£*o 


(39b) 


(39c) 


(39d) 


(40a) 


(40b) 


and the nominal system matrices are given above in equation ( 35 ). It 
should be noted that a certain amount of freedom exists in determining 
the above matrices, so that an uncertainty model obtained for a given 
uncertain system is not unique. It should also be noted that in the 
above uncertainty model development, the scaling terms Spj were 
incorporated into the model at the end so as to reduce the number of 
symbolic terms involved in the determination of the P21, P12. «nd Pi 1 
matrices. 


rational uncertainty modeling problem. This is done by obtaining a 
matrix fraction description of the uncertain system, and representing 
the denominator matrix in a feedback loop so as to remove the inverse. 
Tbe numerator and denominator matrices are that multivariate 
polynomial matrices which can be concate n ated together and modeled 
using the multilinear techniques discussed above. Details of this 
approach will be presented in a subsequent paper* 


6. Conclusions 

This paper has summarized previous results in parametric 
uncertainty modeling, and has presented and demonstrated an 
important extension to these results. The extension consists of a 
framework for modeling multiple parametric uncertainties which can be 
represented in a multilinear functional form, and includes a procedure 
for obtaining a minimal (or near minimal) uncertainty model relative to 
a given state space realization of the uncertain system. As discussed in 
the papa, the multilinear framework can also be used to solve the more 
general rational uncertain parameter case, and provides a mechanism 
for significantly simplifying the computational complexity involved in 
determining an uncertainty model for a given uncertain system. Thus, 
many practical problems of interest can be solved within this 
framework. To demonstrate the results of the paper, an example 
problem was presented which consisted of a multivariable third -order 
uncertain system with four uncertain parameters. A minimal (or near 
minimal) uncertainty model was determined for the given stale space 
realization of this system, and the resulting model had a dimension of 
seven. Although two of the uncertain parameters entered into the given 
model as squared terms and as fractions, they were easily modeled 
within the multilinear framework. 

Further work being addressed in this area includes 
evaluating/refining/general iziing this modeling procedure for a wider 
class of problems, automating the generalized modeling procedure, and 
applying the procedure to practical application problems. 
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5. Extension to Rational Case 


The above procedure for solving the multilinear uncertainty 
modeling problem can in fact also be used to solve the more general 
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